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A simple model for the reheating of the universe after inflation is studied in which an essentially 
inhomogeneous scalar field representing matter is coupled to an essentially homogeneous scalar 
inflaton field. Through this coupling, the potential determining the evolution of the inflaton field is 
made time-dependent. Due to this the frequency of parametric resonance becomes time dependent, 
making the reheating process especially effective. 

All fields including the gravitational field are initially simplified by expanding each in terms of the 
respective homogeneity or inhomogeneity. Employing only the lowest order of this expansion, we 
space-average, and introduce all independent averages as new variables. This leads to a hierarchy of 
equations for the spatial moments of the fields and their derivatives. A small expansion parameter 
permits a truncation to lowest order, yielding a closed system of 5 coupled nonlinear first order 
differential equations. 

For a parabolic potential, the energy densities of the matter and inflaton fields oscillate chaotically 
around each other from the end of inflation until they reach extremely small values. The average 
period of preponderance of one of the two continuously increases in this process. We discuss that 
this may provide one clue to a solution of the coincidence problem. 

For a Mexican hat potential we can easily obtain and understand dynamical symmetry breaking. 

PACS numbers: 95.30.Sf, 98.80.-k, 98.80.Bp, 98.80.Cq 



I. INTRODUCTION 

Inflation has become an almost indispensable ingre- 
dient of cosmology because so far it is the only means 
by which several problems of cosmology, like the hori- 
zon, the flatness or the structure formation problem can 
be solved. Inflation can be achieved by an inflaton field 
mO-El' i- e - a homogeneous field which exerts negative 
pressure and, as the ground state of a quantum field, has 
no particles. In its simplest version it is a scalar field 
which, for simplicity, will also be used here. 

On one hand there must be an end to inflation, and on 
the other, at the end or after inflation there must be for- 
mation and heating and/or reheating of matter which, if 
initially present, has been extremely diluted and cooled 
down by inflation. In order for inflation to come to an 
end, the inflaton field must loose much of its energy while 
for its formation and (re-)heating, matter must acquire 
energy. In principle this can happen independently: Mat- 
ter can extract energy from the gravitational field, and 
energy of the inflaton field can go to the gravitational 
field which effectively amounts to a transfer of energy 
from the inflaton to the matter field. These processes 
can be independent even if they occur at the same time. 
Inflation scenarios employing only an independent infla- 
ton field are of this kind. 

The most convincing inflation scenarios appear to be 
those in which energy is transferred from the inflaton to 
the matter field through direct coupling. In this case it 
was shown that reheating is especially effective due to 
the occurrence of parametric resonance 0, 0, IE S @ ' ^ 



a field description for matter is used. The scenario em- 
ployed in this paper is of this type. We observe chaotic 
behavior which is expected, since there is a close rela- 
tionship between chaos and parametric resonance A 
different approach studying this relationship using lattice 
simulations was considered in flQl . 



II. MODEL 

It will be assumed that the dynamics of the inflaton 
field ip is determined by a potential of the form 

V(<p)=V + ^f + X^ (1) 

which is re-normalizable in a quantum field theory. We 
incorporate matter by introducing a scalar field x, though 
this can only be a crude reflection of reality. An interac- 
tion Lagrangian of the form 



is assumed. A potential V ~ V(tp, x) containing the 
interaction C\ is 

Vfa x)=V +a ( X 2 ~xl) ^ + X^ + ^f. (3) 

In the case that all constants Vq, a, Xc, A and m x are dif- 
ferent from zero, the potential V(ip, x) has, in its de- 
pendence on ip, for given x > Xc a parabola-like shape 
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and for \ < Xc a sombrero-like shape (cf. Fig. 0) . The 
term ~p 2 is implicitly time-dependent, enabling dynam- 
ical symmetry breaking as expected from a Higgs field. 
In addition, setting Vo=0, Xc=0 and m x —0, as a second 
case 
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(4) 



will be considered. 

The model thus introduced resembles hybrid inflation 
models in that it employs two coupled scalar fields. How- 
ever, it differs from them because the second field is not 
an inflaton field but an inhomogcncous field for matter. 

Concerning the space dependence it is assumed that 
the inflaton field starts out as and remains an essentially 
homogeneous scalar field while the matter field starts out 
and remains essentially inhomogeneous, 

ip = ip h (t) + e(p{x,t) , x = £Xh(*) +x(x,t) (5) 

with e <C 1 and <p=x=0, where /= J v f d 3 x/V and V is 
a volume large enough to contain many particles of the 
field x- Note that due to the interaction of the fields, p> 
will get an inhomogcncous contribution even if it starts 
out almost completely homogeneous, while the opposite 
holds for X- 

The purpose of this paper is no study of fluctuations. 
We will, therefore, concentrate on average properties, in- 
homogeneities only being attributed to the field repre- 
sentation of matter. They are necessary for obtaining 
positive pressure but are intended to represent a homo- 
geneous matter distribution on large scales. Thus inho- 
mogeneities occur on a much smaller scale than the one 
of fluctuations used for studying structure formation. As 
a consequence, in an expansion with respect to e we shall 
be satisfied with lowest order results. Therefore, and be- 
cause perturbations of the coordinates would only result 
in second order effects we can use the standard co- 
ordinates of an unperturbed space-time for x and t. It 
should be noted that the following approach is neither 
covariant, nor gauge invariant. However, this does not 
matter because of the just mentioned reasons and the 
focus on physically meaningful quantities. 

The lowest order results which will finally be obtained 
from an expansion with respect to e could easily be 
guessed. Nevertheless, for the reason of clarity they will 
be derived explicitly. 



III. DERIVATION OF THE BASIC EQUATIONS 

Owing to the spatial fluctuations of the fields <p and 
X, fluctuations of the gravitational field will occur. A 
time dependent and spatially fluctuating gravitational 
field can be described by a local scale factor 



a{x, t) = ah(t) + ea{x, t) 



(G) 



and a local Hubble parameter [Tl| 

a(x 7 1) 



H(x, t) 



a(x, t) 



(7) 



(with d=da/dt etc.), which can be decomposed according 
to 



H(x,t) = H h (t) + eh(x,t) 



(8) 



with h=0. Assuming minimal coupling, the fields ip and 
X must satisfy modified Klein-Gordon-equations 



(p + 3Htp — 



X + 3H X 



vy | dv(cp, x ) 

a 2 dp> 

V 2 X | dV{p, X ) 
a 2 dx 



= 0, (9) 



0, (10) 



and H (x, t) must satisfy the Raychaudhuri equation [T^|- 

M 



An 

H +H 2 = -^|(P*+Px+ 3 ?V+ 3 Px) 



A (Pv +Px) 



3a 2 (p v +p x +p v +p x ) 



(11) 



In this the densities p v and p x and the pressures p v and 
p x are given by 



p 2 t (V^ 



2 ' 2a 2 
p 2 (Wp) 2 



x 2 , (V X ) 2 



Vy 



2 6a 2 



v v> Px = ~o 



2 2a 2 ^ Kx ' 
(12) 

X 2 (Vx) 2 



2 6a 



~-Vx, 



where, with proper distribution V — V v + V x of the po- 
tential ©, 



2 4 

V v = Vo-ax 2 c ^- + 



v x = a X 2 tt + 



p 2 m 2 x x 2 



2 4 2 2 

From Equation Q to lowest order in e we get 

cih(t) 



H h (t) 



(13) 



Introducing Eqs. © and © in Eqs. (|9^l- lfTU|) and space- 
averaging Eq. © , to lowest order in e we obtain 



<p h + 3H h (p h + a(x 2 -xi)<Pb + X<pl = 0, (14) 
X + 3H h x 2777 + ocip h x + m 2 X = . (15) 



From Eqs. {BJl with © we get 

Pip = P<pQ{i)+£pv{$,t) , P x = P x o( t )+ e Px(^, t ) i 
P v = P V Q(t)+tp<p{x,t) , p x =p x0 (t)+ep x (x,t) 
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which defines p x , p v , p x and 

(p\ x 2 (Vx) 2 

p v o = -Y+V 9 {yh), P x o = Y+^2^y +v x(^,x) , 

(16) 

<Ph x 2 (Vx) 2 

PfO = -j-VipiVh), P x o = Y~6tfJtj~ Vx<ylfh '^' 

With this and Eq. JBJ), space- averaging Eq. yields 

47T 



H h (t)+H£(t) 



3m 



2-(PvO+PxO+3Pv>o+3p x o) ( 17 ) 



to lowest order in e. Writing 



<p h + 3Hh[t)ifh H 5 = 

for Eq. (|14|l and using the equations p<po+Pipo =l Ph and 
P v o=^\^\i+{dV v {ip\ l ) / diph)^ following from Eqs. {TSJ 
one obtains 



^o + 3i/ h (t)(p y0 +P v o) + ^H^ il ^^h =0. (18) 



Writing Eq. JTSJ as 



X + 3i?hX 



V 2 x , dV x (ip h ,x) 



+ 



0; 



multiplying it with x, space- aver aging it and finally using 



xV 2 x = -Vx • Vx 
which follows from partial integration, we get 



d X 2 , OTT M j£ , V X ■ VX , dV x ((p h ,x) ~ _ n 

-,2 tJ.\ Q ~ X u • 



dt2 ■■ 3ff H(*)x a , ^ . g - 



Combining this with the equations 



— - —^L 4- v * • v * 

Px ° ~ <ft 2 + a 2 (i) 
»X 



-flh(*) 



(Vx) 2 , g^h, x), 
a£(t) %>h 



•<Ph 



It means that to lowest order in e the total energy of 
matter and inflaton field is conserved. It is well known 
that together with Eq. 1)20(1 the equation 



k 



<{t) 



8tt 
3m| 



po 



(22) 



with k=l or or —1 is equivalent to Eq. (|17|l whence we 
can replace the latter with Eq. I|22|) . 

To lowest order in e we are thus left with the set of 
equations l(T3 ^l -l|T5 )l and (J22J- They constitute a self con- 
tained system of four equations for the four unknown 
quantities ah(t), Hh(t), iph(t) and x{x,t) which com- 
pletely decouples from the first order perturbations. In 
the following for simplicity we will set ah — ► a and 

Hh n n 

Equations of a similar kind were solved in [3J, |a| by 

means of a Fourier analysis, invoking quantum field the- 
ory and introducing several approximations. In the fol- 
lowing a quite different approach will be presented. This 
allows for a classical treatment which gets along with mi- 
nor approximations and will enable a simple numerical 
treatment. 

The method chosen here consist in replacing Eq. 1)15(1 . 
the only partial differential equation left in the system, 
by a set of ordinary differential equations for averaged 
quantities. In certain aspects it resembles the momen- 
tum method used for solving the Boltzmann equation, 
the main difference being that in this paper moments are 
built in ordinary space instead of momentum space. Like 
in Boltzmann theory, the quality of the approximation 
can be improved by including moments of ever increas- 
ing order - a typical example is provided by the thir- 
teen, twenty-one or twenty-nine moment approximation 
in Plasma Physics |l3 |. A tremendous advantage over 
the momentum method in Boltzmann theory consists in 
a much better quality of the approximation provided by 
a truncation to low orders. This is enabled by the fact 
that each additional moment entering the equations is 
multiplied by the small quantity 1/a 2 whence indirectly 
a moment of order n couples to the moments of lowest 
order with a factor 1 / a 2n . 

Defining 

x = X , x n = (W n x) 2 , y n = (V"i) 2 (23) 
where V 3 x = VAx etc. one obtains 



and p x o+p x o=X 2 +(Vx) 2 /(3a 2 1 (0) following from iJTfljl we 
get 



p^ + 3H h (t)(p^ + p^)- gMg^l y h = , (19) 

Addition of Eqs. (|18J) and 1)19(1 leads to the equation 

po + 3H h (t)(p a + po) = (20) 



with 



Po — P v o + P x o i 



Po=P v >o+P x o- (21) 



x^ = {\7 n x) = 2(V"i) 2 + 2(\7 n x)(V n x) 
Inserting x=x from Eq. (|15|l in this yields 



2t/„ 



- 3H(V n x)CV n x) + ^r(V"a;)(V™+ 2 a;) 
or 



-{aipi + m x ){V n x)(V n x) 
2jfe - 3Hx^ - 2(aipl + m 2 x )x^ 
+ (2/a 2 ) \div((V n + 1 x)(W n x)) - (V^x) 2 
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The average of the divergence is a surface integral divided 
by V and vanishes for V— >oo whence 



x n = 2y n - 2x n+1 /a 2 - 3Hx n - 2{aip\ + m x )x n . (24) 
Analogously one obtains 



Vn 



-6Hy n - x n+ i/a 2 - (atfQ + m 2 x )x r , 



(25) 



For n=0, 1,2,... Eqs. to S etner with HH and 

(|2*2*|l provide an infinite set of equations coupled together 
in each order n through a coupling term of order n+1. 
The coupling terms are multiplied by the factor 1/a 2 
that, in an expanding universe, is getting smaller and 
smaller with time. In addition, it appears plausible that 
they will be continually diminished by smoothing effects 
of friction. Therefore, it can be expected that truncation 
by neglecting the coupling terms to a low order n will 
yield good approximations. This can be tested numeri- 
cally by comparing results obtained from truncation to 
different orders (see Appendix). It turned out that keep- 
ing the n=0 terms only is already good enough. With 
xq—x 2 and yo=x 2 , the equations obtained this way are 



3H(ph + a (xq 



Xc) V?h 



.to + 3Hx + 2(atp 2 + m 2 )xo 



-250 



2/o 



aipl + m 2 x )x 



a 2 










<Ph 2 


P<pO = 


2 






Pxo = 




a = 


Ha. 



3m 2 



in 



2 

2^0 

2 



4 



0. 


(26) 


o, 


(27) 


o, 


(28) 




(29) 


) 


(30) 




(31) 




(32) 



X 2 



Note that due to the presence of ~x~o — X 2 an d Vo 
the system still contains effects of the inhomogeneity of 
X(x,t), in spite of the neglection of all gradient terms in 
the order n — 0. 

ifh carries out damped oscillations in the neighborhood 
of a minimum of the potential V and plays the role of 
a driver. Eq. I|27|l for xo has properties similar to the 
Mathieu equation, the term 2a<p 2 JEo leading to paramet- 
ric resonance. As a result of this and the nonlinearities 
of the system, chaotic solutions must be expected. 

In the early phase of the universe when the total den- 
sity po=p v o+p X Q is very large, the term k/a 2 in Eq. 129|) 
can be neglected. Calculations for low densities will be 
restricted to the case k=Q. Thus for the purposes of this 
paper Eqs. (|26[1 - H31[) build a closed subset independent 
of the quantity a. With fc=0, 



8tt 
3~' 



X 



Xq 

2 ' 
P 



■in 



4 



III r 



Pv ■ 



PxQ 

4 

77?,p 



Px 



™4 



Ph 

mp 
>V , 



Xc 
TOp 



mpt , (33) 



Xc (34) 



and the definitions £7=$(r), Z=X(r), Eqs. (EHJ dSU for 
the variables ip^, Xq and yo are transformed into the set 
of autonomous first order differential equations 



$(r 
U(t 

Px( T 
X(r 
Z(r 



U. 



= - [3HU + a(X - x 2 c )$ + ^$> 3 ~ 
= -6Hp x + 3aHX<S> 2 + aX$U 



= Z 



-3HZ + 4p x - AaX<5> 2 



with 



H =/3\/P<P+Pxi P v =Vo + 



U 



rXc* 



(35) 
(36) 
(37) 
(38) 
(39) 



(40) 



for the dimensionless variables $, U, p x , X and Z. Once 
the system is solved, a{r) is obtained from 



a exp (/? / Pv (t<) + p x (T<) rfr') . (41) 



IV. 



CONDITIONS ON THE PARAMETERS OF 
THE PROBLEM 



1. As is usual in modern cosmology, we shall assume 
that all quantities are restricted by the corresponding 
Planck quantities. For the (normalized) fields this means 
ip < 1 and x < 1, and for the (normalized) potential V < 1. 

2. Presently the energy densities of the matter and in- 
flaton field are almost equal and very small compared 
to 1. On the assumption that the corresponding fields 
are close to a value where the potential has a minimum, 
the latter must be very small. In case (0J the mini- 
mum is zero. In order to avoid fine tuning we assume 
that this is also true in case In this case, the mini- 
mum is obtained for x=0, ip 2 ~ax 2 /X and has the value 
V min =V -axl/(4\) whence a = 4AF /Xc from V min =0. 

3. A further condition on the parameters of the prob- 
lem originates from the requirement that the observed 
inhomogeneities of the universe evolve from quantum 
fluctuations through inflation. From this, for the poten- 
tial V=Xp 4 /A in [l| the condition A«10 -14 was derived. 
Since matter is getting extremely diluted through infla- 
tion, during this process x~0 whence in case (@J) |x|^C|</?|, 
l/«A<p 4 /4 and thus A~10~ 14 . During inflation x~0 holds 
also for Eq. (J3J) whence V~Vq— ax 2 p 2 /2+Xip 4 /4. In this 

calculation following the lines of that, for rea- 
son of brevity, cannot be presented here yields A«2-10~ 17 
for Vb=1.25-10~ 18 . 



V. INFLATON FIELD MODEL 

The case of the parabola shaped potential (0J is es- 
sentially the scenario of Linde's chaotic inflation pj. 



It is known that for the small value A 
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FIG. 1: (a) Time evolution of p v , p x and 
p during slow roll and first reheating, starting 
with p x =I0" 15 <p ip =0.I5-I0" 11 . Below: Evo- 
lution of the inflaton field $ (b) and mat- 
ter field X (c) over the same time interval. 
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FIG. 2: Time evolution of p v , p x and p in dif- 
ferent time intervals with different time scales, r be- 
ing reset to at the beginning of each: (a) be- 
fore inflation, start with / 9 x =10 _15 <C / o ¥ ,=0.15-10 -11 
(the same as in Fig. QJ, (b) much later, start with 
p x =l(r 30 >p v =lCr 34 , (c) start with / o x =p ¥ ,=5-10 _4x . 



extreme inflationary expansion is obtained if the to- 
tal density p starts at about the Planck density 0. 
Therefore for numerical calculations a much lower ini- 
tial value p(0)«0.16-10 -12 was chosen, leading to an in- 
flationary enhancement of the scale factor a by 35 or- 
ders of magnitude only. In order to obtain inflation, 
P<p0 = — V v ((fh) < must be satisfied which was 

achieved by setting y>h(0) = or U(0) = 0. With this, 
according to Eqs. (f4*U|) 



*(0) 



WO) 

A 



1/4 



For matter, between pressure and density a relation 
Pxo = IPxQ holds where 7 = 1/3 in the highly relativistic 
regime and 7 <C 1/3 in the non relativistic regime. With 

a(p^Xo/2 this 
- 7) and Xq = 



Eq. (|3U|l for m x — and — 

relation leads to = (1 + -f)aipfr^/(l 
(1— 7)pvo/(c*¥'h) or i expressed in the variables introduced 



Eqs. <SH3U, 



X 
Z 



(1-7)^ 



a$ 2 
1 + 7 



1-7 



a$X . 



(42) 
(43) 



These conditions were used as initial conditions for X 
and Z. The exact value given initially to 7 turned out to 
be of no importance in the long run, so we used 7 = 1/3. 

For p x various initial conditions were considered, from 
essentially no matter to equipartition between matter 
and inflaton field energy. The results are discussed in 
the following. 

p x is first dramatically reduced by inflation (Fig.^and 
Fig. [3(a), insert). When $ starts to oscillate about $=0, 
reheating sets in and p x is raised. The frequency of the 
heating oscillations is coupled to the one of the oscilla- 
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tions of 4>(i) and is approximately twice as large. The 
largest p x reached after reheating is a function of a. The 
condition that p x — p v at some point in the evolution 
of the Universe yields the constraint 1.2A < a < 3 A . For 
the numerical calculations a=2 A was used. After p x and 
p v coincide they start to oscillate in a chaotic manner 
around each other. 

A phase portrait of ip shown in Fig. [S] exhibits typi- 
cal chaotic structure. Chaotic behavior arises naturally 
when two or more scalar fields interact with each other 
in an expanding Universe 0, 0, but it is usually 
of a transient type. The observed chaotic oscillations in 
the present model, on the contrary, are persistent. Even 
after quite large excursions p x (t) and p v (t) always find 
back to each other. This kind of intermittent coincidence 
holds true for the whole decay of p=p x +Ptp from the very 
large values before inflation down to very small values like 
those of today. 

It can qualitatively be understood as follows. From 
equations (|35[) - l|37[l and (|40|l one easily derives 

p v (r) = -6H Pv + (3/2)AiJ$ 4 - aX$U (44) 

and 

Px) = -6H(p v -p x )-3aX<P 2 

CIT 

+ ^AiJ$ 4 - 2aX<PU. (45) 

The first term on the right hand side of Eq. (|45() is always 
leading to a decrease of p v — p x . The second term is < 
and the third is > so the two are counteracting, and 
numerically their sum turns out to be smaller than the 
first term. The last term is oscillatory and is dominant 
according to the numerical evaluation. It is a forcing term 
that leads to oscillations of p v — p x because according to 
Eqs. (|37[) and (|44|l one half of it increases p v when the 
other half decreases p x and vice versa. A demonstration 
of this behavior is given in Figs. OH which show results for 
P<p( t ) - Px (t) and 2oX(t)$(t)U(t). 

Some insight into the oscillations involved in the pro- 
cess is obtained by the following consideration. For 
Xc = combining Eqs. and yields 

$(r) = -3i?$ - al$ - A$ 3 . (46) 

This is an equation for the oscillations of <I> in which 
— (aXQ + A<I> 3 ) is the driving force and — 3H<b is a fric- 
tion term. After a few oscillations through $ = the 
dominant term becomes -aX$, and from 

i>(r) w -aA$ 

the oscillation frequency 

lj$ « JaX (47) 

is obtained. From Eqs. l|3"£ ]) -l|3l? )l we get 

X{t) = -3HX + Ap x - 4aX$ 2 . 




FIG. 3: Density difference /%>(t) — Pxi T ) an d driving 
term 2aA(r)$(r)t/(r) of equation (|43|l . The same set- 
tings and initial conditions as in figure 121(c) were used, 
but only about half of the time-interval was covered. 

Again this is an equation for oscillations, this time of 
A(t), with — AaXtfi 2 being the driving term, and an ap- 
proximate expression for the frequency of the oscillations 
is 

oj x ~2\$\V^- (48) 
According to the numerical calculations, on the average 
ux_ 2J£| 

w$ Va 

is smaller than 1 whence the oscillations of X are slower 
than those of <&. In consequence, in the short run the 
frequency w$ « of the <!> oscillations varies between 
smaller and larger values as can be clearly seen in Fig- 
ure 0] and is also reflected in Fig. |3J In the long run X 
decreases like all other variables of the system of Eqs. 
(13511 14UI) . Therefore, as time goes on, the average fre- 
quency of the $ oscillations and through this also that of 
the oscillations of p v and p x around each other is slowly 
decreasing. 

Obviously the chaotic nature of the system is still quite 
structured, because roughly the two frequencies lox and 
w$ remain visible as organizing factors. It is well known 
that in linear systems with two incommensurable fre- 
quencies ergodic phase portraits are obtained. It appears 
that the chaoticity of the present system is caused by a 
permanent change of the ratio uq/ujx- In this process, 
due to the overwhelming abundance of irrational num- 
bers, the change between different measures of incom- 
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FIG. 4: Approximate frequencies w$ » vol 
and wx ~ 2 1$| y^a, for the same settings and 
initial conditions as used in figure 121(c), but 
only about half of the time-interval was covered. 



mensurability leads to a certain amount of unpredictabil- 
ity. 

Because it was not possible to calculate the full evolu- 
tion from inflation to the present state, the time-averaged 
coincidence of p v and p x claimed above was tested piece- 
wise in more than ten different p-regimes from very large 
to very small, starting with many different values of the 
ratio Pxl Pv m each. Fig. El shows 3 typical examples. 

For the evolution of g after inflation the approximate 
validity of ga v =const was found with i/as4 in (a), i/«3.9 
in (b) and ^«4 in (c) etc., i.e. the usual relativistic equa- 
tion of state, pa 4 =const, is a good approximation to the 
evolution of the (rather smooth) total density g. 

VI. REMARKS ON THE COINCIDENCE 
PROBLEM 

In Ref. fl 6 | an idea for a solution of the coincidence 
problem [lj| was elaborated: The inflaton field energy, 
called dark energy there, has periodically dominated in 
the past, giving a finite probability to its observed pre- 
ponderance today. The oscillations are achieved by sinu- 
soidal modulations of an exponentially decaying potential 
for the scalar field representing dark energy, i.e. by the 
assumption of a periodic forcing. 

Our model leads to a similar behavior, only that the 
oscillations are not pur ely periodic but slightly chaotic. 
In contrast to Ref. |16| they are not implied by an im- 




FIG. 5: Phase portraits of the inflaton-field <j>: (a) 
the same time interval and same initial settings as 
in Fig. n were used, only the slow roll phase is 
omitted; (b) continues where (a) stopped and cov- 
ers the rest of the time interval used in Fig. [21 (a) - 
the chaotic nature of the oscillations is easy to see. 



posed forcing but caused by the oscillations of the infla- 
ton field $ around $ = and its interaction with the 
matter field \- The mechanism is the following: The 
oscillations of $(i) cause heating oscillations by which 
again and again the energy density of matter is temporar- 
ily raised above the energy density of the inflaton field. 
Since with decreasing amplitude of the <& oscillations also 
the frequency is getting smaller and smaller, the periods 
with a preponderance of either matter or inflaton field 
become increasingly longer. 

In our model the alternating preponderance of matter 
and inflaton field was observed over a density range from 
10~ 12 to 10 -115 with the same set of parameters a and 
A. However, in spite of its continuous increase at low 
densities the average period of preponderance is still too 
small by orders of mag nitude. With A = 10" 120 instead 
of A = 10~ 14 a reasonable period (small fraction of the 
life time of the universe) can be achieved, but with this 
small value the behavior at large densities is no longer the 
needed reheating scenario, which was found only down to 
values A> 10~ 40 . At first glance it would appear that a 
potential V($) that decays like 10~ 14 $ 4 at large $ and 
like 10 -120< f> 4 at small <E> would lead to a solution. It 
turns out, however, that an effective reheating at large 
densities requires a large value of A even around <f> = 0. 
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FIG. 6: Potential © for 4 different values of Tqx that de- 
crease in the shape sequence 1,2,3,4. During pre-inflation 
shapes 1 and 2 are active, during inflation shape 3. 4 
is the final shape when matter is already very diluted. 
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FIG. 7: Time evolution of p v and p x for 

the potential © that evolves according to Fig. |2 



At second glance a dynamic A = A(i) would appear 
to do it, and X(t) ~ H 2 (t) would look like a good 
choice, coupling the potential directly to the energy of 
the gravitational field. However, in a reasonable dynam- 
ical evolution pressure and density should be connected 
by p^j — 7/5^o and in consequence, Eq. (|42l) should hold. 
From this and l|47[) we get 

~ Vfa ~ H > H 

W * ~ |$| ~ |$| - max |*f 

While at later times ui$ w i? would have the desired 
order of magnitude because is approximately the 

age of the universe, the small factor 1/|$| spoils any 
attempt to obtain reasonably small frequencies. 

From our studies it can be concluded, that the oscilla- 
tions of the inflaton field and its interaction with matter 
can keep the energy densities of the two close together 
over the full range of densities from inflation to present 
day densities, with alternating preponderance of one of 
the two. Reasonable times for the duration of the pre- 
ponderance at present day densities cannot be obtained 
from our model. So oscillations of the inflaton field and 
periodic reheating of matter through a coupling between 
matter and inflaton field may be one clue to the solu- 
tion of the coincidence problem but cannot provide a full 
solution of it. The simple model for describing matter 
adopted in this paper is certainly not appropriate for 
later stages in the evolution of the universe, and a better 
adapted model might improve the situation. 

VII. HIGGS FIELD MODEL 

For the potential calculations were done 

with F =1.25-10- 18 , A=2-1CT 17 , a=A/10, Xc=V5, 
m^,=l-10~ 19 and with the initial values p x —7.5-l0^ 16 , 
/5<p=3T0 of the densities. First the potential assumes 
shape 1 of Fig. and ip starts at a value for which V 
is well above its minimum at V=0. As ip rolls slowly 
towards this, p x decreases whence V gradually assumes 



a flatter shape like 2. When ip reaches the region around 
ip=0, p x may be enhanced through oscillations of ip in 
the potential well but is then again diluted by ongoing 
inflation. Finally xq=x 2 becomes smaller than \\ and 
V(p) assumes a sombrero-like shape as 3. Then ip starts 
to roll downhill towards a new minimum appearing at 
a ip^O- Oscillations around this lead to an increase of 
p x . Due to the change of shape of the potential the 
oscillation frequency changes and thus parametric reso- 
nance occurs at many different frequencies, leading to an 
especially effective reheating. Asymptotically, shape 4 
will be assumed and ip will settle at the minimum V=0 
obtained from Eq. (|3J) for x=0. The dynamics can be 
separated into four qualitatively different steps (Fig. [7J . 
1. Pre-inflation: Slow roll in a potential of shape 1-2 
towards tp=0 with possible reheating through oscillations 
around ip=0. 2. Inflation: Slow roll from ip=0 in a 
potential of shape 3 towards a new minimum at <p=/=0. 
3. Reheating: Oscillations around the new minimum. 4. 
Friedmann-like evolution after the settlement of p at the 
new minimum. These four steps are found for a wide 
variety of parameters and initial conditions. 

A model combining the properties of the two cases con- 
sidered, chaotic oscillations with alternating dominance 
of matter and inflaton field as well as dynamic symmetry 
breaking, should be obtainable by employing two infla- 
ton fields of type ip = ph(t) + stp(x, t), one moving in a 
parabola- and the other in a sombrero-like potential, and 
both being coupled to the matter field. 

VIII. CONCLUSIONS 

The main result of this paper consists in a set of equa- 
tions for the interaction between a scalar inflaton field, a 
scalar matter field and the gravitational field, obtained by 
averaging processes. Due to the averaging no specific as- 
sumptions about the inhomogeneities of the matter field 
like specifying any initial conditions must be made. The 
set of equations can be applied to pre-inflation, inflation 
and the reheating after inflation. It contains the heating 
effects of parametric resonance and is nevertheless rather 
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FIG. 8: Time evolution of p v , p x and p computed to or- 
der n = 1. The same time interval and initial conditions 
as in Fig. |21 (a) were used, except for a and xT which was 
set to di = nip 1 and y^af/lO respectively. Comparison 
with Fig[21 (a) shows a similar evolution of the densities. 

easy to handle. However, it should be noted that the 
presented approach is neither covariant, nor gauge in- 
variant. Therefore one has to take great care to compute 
only physically meaningful quantities, as in the studied 
cases. 

Processes like the dynamic symmetry breaking of a 
Higgs field, the theoretical description of which usually 
requires quite some effort, arc easily obtained and under- 
stood. An attempt to gain some insight into the coin- 
cidence problem was made. Repeated periods of matter 
reheating can keep the energy densities of the inflaton 
and the matter field close together for all times until to- 
day and lead to an alternating dominance of one of them. 
The preponderance of the inflaton field observed today 
would thus be a question of chance and be preceded by 
a preponderance of matter. However, the specific model 
employed for describing this process did not yield reason- 
able durations of preponderance over the full evolution 
from inflation until today. We therefore conclude that it 
may provide one clue to the understanding of the prob- 
lem while other clues may still be missing. 
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IX. APPENDIX 

The quality of the approximation provided by equa- 
tions l|35[) - (|4U[l can easily be tested by truncation to 
higher orders. Truncating to order n = 1 enhances our 
set l|55|) - (PUl) by two equations and the two variables X\ 
and yi as 

Xq = 2yo — 2x~[/ a 2 — SHxfi — 2(a(p h + m 2 ,)^, 
Vo = -6Hlfi -xT/a 2 - (a(pl + m 2 ,)^ ; 
5TT = 2yT — ?>H~x~i — 2(aip 2 1 + ttl^Jxi , 
m = -6i?yT - (axph + m^jTi . 

In addition, one needs to compute the time evolution of 
the scale factor a by means of equation (|32|) simultane- 
ously, since it enters the equations now. 

In view of their chaotic character a detailed comparison 
of the fine structure of the solutions obtained with trun- 
cation to order n = and to order n = 1 is not meaning- 
ful. However, the average behavior in the long run must 
be comparable, and indeed all computations yielded the 
same qualitative and essentially the same quantitative 
behavior - compare Fig. [5] (a) with Fig.[H]as an example. 
No crucial dependence on the initial values of aTj" and 
yi was observed, so we choose Xq = (1 — j)p x /(a$ 2 ), 

yE = (1 + j)p x /2,x~[ = (1 + 7)Px/ 2 and 2/i = °> witn 
7 = 1/3 in Fig. |HJ As an initial value for the scale fac- 
tor we choose a = trip 1 , corresponding to a Planck sized 
patch of the Universe. 

By the same token the Higgs field model of section IVTTI 
was checked up to second order and no essential difference 
of the long run behavior was found either. 

From this comparison it can be concluded that trun- 
cation to the order n = already provides a very good 
approximation. 
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